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Abstract 

The experimentally measured elliptic (V2) and hexadecapole (V4) flow of charged particles as a 
function of transverse momentum (pj) at midrapidity in Pb-Pb collisions at -\fsmi = 2.76 TeV 
are compared with the relativistic viscous hydrodynamic model simulations. The simulations are 
carried out for two different initial energy density profiles obtained from (i) Glauber model, and 
(ii) Color Glass Condensate (CGC) model. Comparison to experimental data for 10-20% to 40- 
50% centrality, shows that a centrality dependent shear viscosity to entropy density (rj/s) ratio 
with values ranging between 0.0 to 0. 12 are needed to explain the v% data for simulations with the 
Glauber based initial condition. Whereas for the CGC based initial conditions a slightly higher 
value of rj/s is preferred, around 0.08 to 0.16. From the comparison of the V4 simulated results 
to the corresponding experimental measurements we observe that for the centralities 20-30% to 
40-50% the rj/s values lies between 0.0 to 0.12 for both the initial conditions studied. The tj/s 
values obtained from our studies for Pb-Pb collisions at -\fstw = 2.76 TeV are compared to other 
studies which uses both transport and hydrodynamic approaches. 



1. Introduction 

One of the interesting aspects of the high energy heavy-ion collisions at the Large Hadron 
Collider (LHC) is that, a deconfined state of matter consisting of quarks and gluons, if produced 
in the collisions, exists for a sufficiently long duration so as to enable the study of its transport 
properties. The shear viscosity (77) is one of the transport coefficients that can be estimated 
by comparing the experimental data to the theoretical calculations based on either transport 
or hydrodynamic formulations. Particularly the experimental measurements of the azimuthal 
anisotropy of produced particles are compared to corresponding theoretical simulations to obtain 
an estimate of shear viscosity to entropy density ratio (rj/s). Where the azimuthal anisotropy 
measurements in the experiments are characterized in terms of various coefficients of the Fourier 
expansion of the azimuthal distributions of produced particles relative to the reaction plane. The 
second coefficient, v'2, is called as the elliptic flow and is the most widely used observable to 
extract the rj/s value. Such approaches have been successfully applied to the experimental data 
for the Relativistic Heavy-ion collider (RHIC) energies. The extracted rj/s values at RHIC are 
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found to be lie between 1 - 5 x l/47r|l],0000|6l0]. n/s = 1/4tt is the lower bound predicted 
by P. Kovtun et al. [8] (known as KSS limit) for strongly coupled quantum fluids using the string 
theoretical approaches. 



At the LHC energies for Pb-Pb collision at y/s^N — 2.76 TeV one expects the initial tem- 
perature attained in the collisions to be higher compared to that at RHIC for Au-Au collisions at 
y*AW = 200 GeV. Hence it would be interesting to find out the value of rj/s for the QCD matter 
produced at a higher temperature at LHC compared to RHIC. The estimates so far carried out 
based on both transport theory |0 [Tel and viscous hydrodynamic simulations ifTH [T2I [l3l [l4Tl 
have shown that the values of rj/s at LHC energies are similar to those obtained at RHIC. The 
hydrodynamic approaches used included 2+1D simulations with smooth Glauber based initial 
conditions and 3+ ID simulations with fluctuating initial conditions. These approaches mostly 
concentrated on explaining the measured v'2 for certain collision centrality classes. One of the 
dominant source of uncertainties in such estimates is the lack of complete knowledge of the ini- 
tial conditions in the heavy-ion collisions. In this paper, we discuss results from a 2+1D smooth 
viscous hydrodynamic simulations with two different initial conditions, one based on Glauber 
model and other based on CGC approach. These simulations are then compared to centrality 
dependence of experimentally measured V2 as well as new data on hexadecapole flow (V4) as 
a function of transverse momentum from ATLAS lfl5ll and ALICE experiments lfl7ll at LHC. 
Our estimated values of rj/s based on simulations with two different initial conditions and two 
observables, v>2 and V4 are then compared to other existing results in the literature. 

The paper is organized as follows. In the next section we discuss the details of our 2+ ID 
viscous hydrodynamic simulations. These includes the initial conditions, an equation of state 
based on the lattice QCD and hadron resonance gas, and the freeze-out criteria. In section 3 
we present the comparison of our simulation results to the corresponding experimental data at 
midrapidity for Pb-Pb collisions at y/s/w = 2.76 TeV. These includes the centrality dependence 
of the charged particle pseudorapidity density (dNchldrj), V2(pr), and v^pt). In section 4 we 
present the predictions for charged hadron pj spectra from our simulations for various collision 
centralities in Pb-Pb collisions at y*MV = 2.76 TeV. In section 5 we compare our estimated rj/s 
values to other calculations. Finally we summarize our work and present a brief outlook in the 
section 6. 

2. 2+1D viscous hydrodynamic simulations 

The hydrodynamic simulation results presented here are based on the second order causal 
Israel-Stewart formalism for the evolution of a viscous fluid using the 2+1D viscous hydrody- 
namic code "AZHYDRO-KOLKATA" Hd. Shear viscosity is the only dissipative process 
considered in our present study. We assume a net-baryon free plasma is formed in Pb-Pb colli- 
sions at midrapidity at y^jv - 2.76 TeV, which is supported by the value of anti-proton to proton 
ratio being close to unity [20]. The code solves the following energy-momentum conservation 
equation along with the relaxation equation for shear viscosity. 
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Where T MV — (s + p)u^u v - pg^ v + n^' is the energy-momentum tensor, s, p, u, and n^ v aie the 
energy density, pressure, fluid velocity and shear tensor respectively. D = u^d^ is the convective 
time derivative, V <m m v> = \(Wu v + W) - \(d p u fl )(g ,lv - u^u y ) is a symmetric traceless tensor. 
rj is the shear viscosity and t„ is the corresponding relaxation time. Assuming boost-invariance, 
the equations are solved in (t, x, y,rj) coordinate system. 

2.1. Initial conditions 

The initial conditions used here includes the initial energy density profile in the transverse 
plane (e(x,y)). The e(x,y) is obtained using two different models - Glauber and CGC, the initial 
time (t,), the transverse velocity profile (v x (x, y), V y (x, y)), shear stresses in the transverse plane 
{n^ix, y)) at t,. The r, value is taken as 0.6 fm. The temperature independent rj/s values which 
are used as an input to the viscous hydrodynamics simulation are considered to be 0, 0.08, 0.12, 
and 0.16. 

The two component Glauber model based initial condition for the transverse energy density 
is constructed as follows. At an impact parameter b for Pb-Pb collisions at yi^w = 

2.76 TeV, 

the transverse energy density is obtained as, 

e(b,x,y) = e [(l - x h )^^-(b,x,y) + x h N coU (b,x,y)] (3) 

where N par ,(b, x,y) and N co u(b, x,y) are the transverse profile of participant numbers and binary 
collision numbers respectively. The eo corresponds to the central energy density for b — col- 
lision and does not depend on the impact parameter of the collision. The parameter Xf, is the 
hard scattering fraction. Both eo and Xf, are fixed to reproduce the experimental charged hadron 



multiplicity density at midrapidity {dN c hldrj) of 1601 ± 60 for the central 0-5% collisions [21] 
for each value of rj/s. No further adjustment of the parameters are done for other collision cen- 
tralities. The values of eo for different values of rj/s are shown in TableQ] The N part (b, x,y) and 



N co n(b, x, y) values are obtained using an optical Glauber model calculation 12211 . The inelastic 
nucleon-nucleon cross section is taken as 72 mb. 



Table 1: Values of eo used in Glauber, CGC model and the normalization constant C used in CGC model for initial 
transverse energy density. 



rj/s 


e (GeV//m j ), Glauber 


C (GeV//m 1/3 ) CGC 


e (GeV//m J ), CGC 


0.0 


112 


0.033 


77.7 


0.08 


98 


0.030 


70.6 


0.12 


88 


0.027 


63.5 


0.16 


78 


0.024 


56.5 



The CGC based initial transverse energy density profile is obtained as follows. We follow 
references J23 , 24 ] and consider that the initial energy density can be obtained from the gluon 
number density through the thermodynamic relation, 



e(Ti,x T ,b) -Cx 



dN, 
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where J^* dY is the gluon number density evaluated at central rapidity Y — and the overall 
normalization C is a free parameter. C is fixed to reproduce the experimentally measured charged 
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particle multiplicity density at midrapidity. The values of C used in the simulations for different 
input values of 77/ s and the corresponding central energy densities are given in Table Q] The 
details of the implementation of this approach can be found in Refs ffl, |25ll26ll . 

Finally the shear stresses are initialized to their corresponding Navier-Stokes estimates 
for the boost invariant velocity profile, n xx - n yy = 2t7/3t,-, and = 112711 . We have used 
T n - 3Tj/4p (where 77 and p are the shear viscous coefficient and pressure) in our simulation, 
which corresponds to the kinetic theory estimates of relaxation time for shear viscous stress for 
a relativistic Boltzmann gas M28I1 . The initial values of v x (x,y) and v y (x,y) are taken to be zero. 

2.2. Equation of state 

The equation of state used here is a combination of high temperature phase from the Lattice 



QCD result with cross-over transition at temperature T c — 175 MeV [29] and the low tempera- 
ture phase is modeled by hadronic resonance gas, containing all the resonances with M res <2.5 
GeV 113011 . Entropy density (s) of the two phases are joined at T — T c = 175 MeV by a smooth 
step like function. The thermodynamic variables pressure (p), energy density (e), etc. are then 
calculated by using the standard thermodynamic relations. 

p(T) = f s(T')dT' (5) 



s(T) = Ts(T)-p(T). (6) 

2.3. Freeze-out condition 

We use the Cooper-Frey algorithm at the freezeout to calculate invariant yields of the hadrons 113 ill 
. The freezeout temperature which is a free parameter in the hydrodynamics simulation is taken 
as 77= 130 MeV. 

The freezeout distribution function in presence of viscous effects gets modified and can be 



approximated as tl8 ] 



fneq(x,p) = f eq (x, p)[l + <t>(x,p)], (7) 

where <p(x, p) « 1 is the corresponding deviation from the equilibrium distribution function 
feq(x, p). The non-equilibrium correction <p{x, p) can be calculated in Grad's 14 moment method. 
According to which the viscous correction terms up to a quadratic function of the four momentum 



p M 02ll33ll . As only shear stresses are considered, cf>(x, p) has the following form 

<l>{x,p) = e iiv p»p v , (8) 

where 

s " v = 2( e + p)r- n ^ (9) 

The correction factor thus increases with increasing values of shear stress at freezeout. The 
correction term also depends on the particle momentum. The Cooper-Frey formula l3~ltl for a 
non equilibrium system is 



in where g is the degeneracy of the particle considered and dL^ is the normal to the elemental 
112 freeze-out hypersurface. 
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Table 2: Comparison of experimentally measured dN c i,/dr] at midrapidity for Pb-Pb collisions at y*NN = 2.76 TeV for 
various collision centrality |21] to those obtained from a 2+ ID viscous hydrodynamic simulations with Glauber and 
CGC based initial conditions. 



% cross section 


Experiment 




dN c h/dr] 
Glauber 










77/s = 0.0 


77/s = 0.08 


77/s = 0.12 


77/s = 0.16 


0-5 


1601 ±60 


1569 


1572 


1567 


1565 


10-20 


966 ± 37 


959 


968 


966 


967 


20-30 


649 ± 23 


637 


647 


649 


651 


30-40 


426 ± 15 


414 


424 


427 


428 


40-50 


261+9 


248 


256 


260 


260 


% cross section 


Experiment 




CGC 










77/s = 0.0 


77/s = 0.08 


77/s = 0.12 


77/s = 0.16 


0-5 


1601 ±60 


1577 


1629 


1623 


1626 


10-20 


966 ± 37 


974 


981 


985 


988 


20-30 


649 ± 23 


640 


666 


668 


670 


30-40 


426 ± 15 


424 


443 


446 


447 


40-50 


261 ±9 


259 


274 


276 


276 



113 3. Comparison to experimental data 

114 The experimental data used for comparison to the simulated results are from the ALICE and 

115 ATLAS experiments at LHC. The observables are charged hadron multiplicity density {dN c i,/dri) 

116 at midrapidity ll2lll . elliptic flow lfl5ll . and hexadecapole flow lfl7ll as a function of pj measured 

117 at midrapidity for Pb-Pb collisions at y/s^N = 2.76 TeV. All results are for collision centrality 
us 10-20%, 20-30%, 30-40%, and 40-50%. For charged hadron yields we also present results for 

119 0-5% collision centrality. This centrality is used to fix the initial energy density parameter in our 

120 simulations. 

121 The V2(pr) data is from the ATLAS experiment lfl5tl and are obtained using the formula 

122 V2 = (cos(2((f> - ¥2))) after correction of the event plane resolution. Where <p is the azimuthal 

123 angle of the charged hadrons measured within | 77 | < 1 and ¥2 is the second order event plane 

124 constructed using a forward calorimeter detector in 3.2 < 77 < 4.9. The rapidity gap between the 

125 detectors used to measure the V2 and ¥2 ensures absence of significant A77 dependent non-flow 

126 correlations, which are also absent in our hydrodynamic simulations. T*2 is zero, i.e, the impact 

127 parameter vector 'b' is along x axis. 

128 The V4(pr) data is from the ALICE experiment [ 17] and are obtained using the formula V4 = 

129 (cos(4(4> - v r , 4))) after correction of the event plane resolution. Where (f> is the azimuthal angle of 

130 the charged hadrons measured within 1 77 1 < 0.8 and ¥4 is the fourth order event plane constructed 

131 using two scintillator arrays of detectors covering -3.7 < 77 < —1.7 and 2.8 < 77 < 5.1. Like 

132 V2(pr), the rapidity gap between the detectors used to measure the V4 and ¥4 ensures absence 

133 of significant A77 dependent non-flow correlations, which are also absent in our hydrodynamic 

134 simulations. ¥4 is zero for our simulation and impact parameter 'b' is along x axis. 
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Figure 1: (Color online) Elliptic flow of charged hadrons as a function of transverse momentum at midrapidity for 
Pb-Pb collisions at yj'NN = 2.76 TeV. The open circles corresponds to experimental data measured by the ATLAS 
collaboration 1 15] . The lines represent results from a 2+1D relativistic viscous hydrodynamic model with a Glauber 
based (upper panel) and CGC based (lower panel) initial transverse energy density profile and different rj/s values. 
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135 3.1. Charged particle pseudo rapidity density 

136 The initial energy density profiles are fixed in the simulations by matching the experimen- 

137 tal measured dN c h/dj] for only 0-5% central Pb-Pb collisions at midrapidity for V^nn = 2.76 

138 TeV I21I1 . No adjustment of the parameters are done to match the dN c h/dT] for the rest of the 

139 collision centralities 10-20%, 20-30%, 30-40%, and 40-50% studied in this paper. The table El 
mo shows the comparison between dN c i,/dri at midrapidity as measured by ALICE for Pb-Pb colli- 

141 sions at y^NN = 2.76 TeV to our simulation results for various values of rj/s for both Glauber 

142 and CGC based initial conditions. We observe for the collision centrality classes studied, there is 

143 good agreement between the experimental measurements and the simulated results for both the 

144 initial conditions. 

145 3.2. Elliptic flow 

146 Upper panel of Fig Q] shows the V2 as a function of pr for charged hadrons at midrapidity (| 

147 77 1< 1) for Pb-Pb collisions at ysNN = 2.76 TeV. The results are shown for four different collision 
we centralities (10-20%, 20-30%, 30-40%, and 40-50%). The open circles are the experimental data 

149 from the ATLAS collaboration [15]. The simulated results are from the 2+lDrelativistic viscous 

150 hydrodynamic model with a Glauber based initial transverse energy density profile. The black 

151 solid, red dashed, blue dotted, and magenta dot-dashed lines corresponds to calculations with rj/ s 

152 = 0.0, 0.08, 0.12, and 0.16 respectively. The experimental data for 10-20% collision centrality 

153 is best described by ideal fluid (rj/s = 0.0) simulation results for most of the pj range for which 

154 comparison has been done, the experimental data for 40-50% collision centrality prefer a tj/s 

155 between 0.08 and 0.12. 

156 Lower panel of Fig. [1] shows the same results as in the upper panel but now the simulated 

157 results corresponds to 2+ ID viscous hydrodynamic calculations with a CGC based initial trans- 

158 verse energy density profile. Comparison of the simulated results in the upper panel and lower 

159 panel for the same collision centrality shows that the viipr) for CGC based initial condition is 

160 larger compared to corresponding results from Glauber based initial condition. This can be un- 

161 derstood from the fact that CGC based initial condition leads to a higher value of momentum 

162 anisotropy compared to Glauber based initial condition [ l].We find from the comparison of ex- 

163 perimental data to simulations based on CGC initial conditions that the v^ipr) data for all the 

164 collision centralities studied are best explained for with rj/s between 0.12-0.16. In contrast to 

165 Glauber based initial conditions we find the CGC based initial conditions prefer a slightly higher 

166 value of 77/. s for the LHC data. 

167 3.3. Hexadecapole flow 

168 Fig. El shows the hexadecapole flow (V4) as a function of transverse momentum (pr) for 

169 charged hadrons at midrapidity in Pb-Pb collisions at V^nn = 2.76 TeV. The results are shown for 

170 four different collision centralities (10-20%, 20-30%, 30-40%, and 40-50%). The open circles are 

171 the experimental data from the ALICE collaboration lfl7ll . The upper panel shows the simulated 

172 results for ideal and viscous fluid evolution using Glauber based initial condition, while the 

173 simulated results for the CGC based initial conditions are shown in the lower panel. We find that 

174 \'a{pt) from ideal hydrodynamic simulations with Glauber based initial conditions is smaller 

175 than the experimental data for 10-20% collision centrality and agrees with the experimental data 

176 at lower pj for 20-30% collision centrality. For the more peripheral collisions of 30-40% and 

177 40-50% centrality, the comparison of the simulated results to the experimental data shows that 

178 the preferred rj/s lies between 0.0 and 0.08. A similar conclusion is also obtained from the 

179 comparison between simulated results with CGC based initial condition and experimental data. 
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Figure 2: (Color online) Hexadecapole flow of charged hadrons as a function of transverse momentum at midrapidity 
for Pb-Pb collisions at y.«NN = 2.76 TeV. The open circles corresponds to experimental data measured by the ALICE 
collaboration 1 17]. The lines represent results from a 2+ ID relativistic viscous hydrodynamic model with both Glauber 
based (upper panel) and CGC based (lower panel) initial transverse energy density profile for different r]/s values. 



8 




1 2 3 4 1 2 3 4 

Transverse momentum, p T (GeV) 

Figure 3: (Color online) Invariant yield of charged hadrons as a function of transverse momentum at midrapidity for 
Pb-Pb collisions at yj'NN = 2.76 TeV. The open circles corresponds to experimental data for 0-5% centrality collision 
measured by the ALICE collaboration I19I . The lines represent results from a 2+ ID relativistic viscous hydrodynamic 
model with a CGC based (left panel) and Glauber based (right panel) initial transverse energy density profile for different 
r\js values. 



so 4. Simulated charged hadron transverse momentum spectra 

181 We have compared our simulated results for charged hadron transverse momentum spectra 

182 with the available experimental data for 0-5% centrality Pb-Pb collision at LHC |l_6]. The com- 

183 parison is shown in Fig[3]. The left panel of Fig[3]is for CGC initial condition whereas the right 

184 panel of the same figure shows the result obtained for Glauber initial condition. We have fixed 

185 the initial energy density of the fluid and the freeze-out temperature by reproducing the expert- 
ise mental transverse momentum spectra of charged hadron for 0-5% collision centrality. So it will 

187 be illogical to extract viscosity by comparing simulation results to the experimentally measured 

188 transverse momentum spectra for 0-5% centrality. 

189 Fig [4] shows the simulated invariant yield of charged hadrons as a function of transverse 

190 momentum at midrapidity for Pb-Pb collisions at -y^NN = 2.76 TeV for four different collision 

191 centralities (10-20%, 20-30%, 30-40%, and 40-50%). The open symbols corresponds to the 

192 simulated results from the 2+1D relativistic viscous hydrodynamic model with a Glauber based 

193 initial transverse energy density profile while the lines corresponds to the calculations with a 

194 CGC based initial condition. The black open circle (black solid line), red open square (red 

195 dashed), blue open triangle (blue dotted) and magenta diamond (magenta dash-dotted) symbols 

196 corresponds to simulations with rj/s = 0.0, 0.08, 0.12, and 0.16 respectively. As we go from 10- 

197 20% collision centrality to 40-50% collision centrality for CGC based initial condition, spectra 

198 become more flatter than for Glauber initial condition. This means that the average transverse 

199 velocity at the freeze-out which determines the slope of the pj spectra is slightly larger for the 

200 fluid evolution with Glauber than those corresponding to the CGC based initial conditions. 
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Figure 4: (Color online) Invariant yield of charged hadrons as a function of transverse momentum at midrapidity for 
Pb-Pb collisions at ysNN = 2-76 TeV from a 2+1D relativistic viscous hydrodynamic model. The open symbols are 
for results obtained using Glauber based initial transverse energy density profile and different r\js values. The lines 
corresponds to calculations with CGC based initial transverse energy density profile. The results are shown for four 
different collision centralities 10-20%, 20-30%, 30-40%, and 40-50%. 
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Figure 5: (Color online) Estimates of T]/s using both transport and hydrodynamic approaches by different groups. The 
results are expressed in terms of multiples of the KSS limit. 
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5. Comparison of our estimate of i]/s with other calculations 

Fig [5] shows the comparison of the estimated rj/s from the experimental data at LHC ener- 
gies from various approaches to one obtained in the current work. Most of the results seems to 
indicate the value of rj/s < 4 times the KSS limit. Two of the approaches shown in the figure 
are based on microscopic transport theory. The remaining are based on the macroscopic hydro- 
dynamic approaches. J. Xu et al. |@], have estimated the specific shear viscosity using A Multi 
Phase Transport model (AMPT). They have explained the experimental data on charged parti- 
cle pseudorapidity density per participating nucleon pair, V2(pr) hi 40-50% centrality, and the 
centrality dependence of V2 for Pb-Pb collision at = 2.76 TeV. S. Plumari et al. Ill Oil , have 

used a parton cascade based approach with a temperature dependent J]/s(T) to simulate V2 at 
LHC energies. They compare the simulated results to experimental measurements of V2(pr) for 
20-30% collision centrality Pb-Pb collisions at -\/Jnn = 2.76 TeV. P. Bozek 111 111 has estimated 
the value of rj/s by using a 2+ ID viscous hydrodynamic model in which a Glauber based entropy 
density initialization was considered. In addition to shear viscosity a bulk viscosity to entropy 
density of 0.04 was also considered in the simulation. Further the freeze-out and resonance de- 
cay contributions were based on THERMINATOR event generator. M. Luzum lfl2ll uses a 2+ ID 
viscous hydrodynamic simulations to carry out shear viscous evolution starting from a smooth 
initial condition for LHC energies. B. Schenke et al. lfl3ll . uses a 3+1D viscous hydrodynamics 
with fluctuating initial conditions to explain the experimental data on V2(pr) and pj integrated V2 
for various collision centralities to estimate the value of r//s at LHC energies. H. Song IU4I1 has 
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221 used a hybrid approach to explain the measured pj integrated V2 as well as pj differential V2 for 

222 various collision centrality in Pb-Pb collisions at -^sjyiv = 2.76 TeV. The hybrid approach com- 

223 bines the macroscopic viscous hydrodynamic description for the QGP fluid with the microscopic 

224 hadron cascade model for the subsequent evolution of the hadronic stage. 

225 6. SUMMARY 

226 In summary, within the framework of Israel-Stewart's 2nd order hydrodynamics, we have 

227 simulated ysAw=2.76 TeV Pb-Pb collisions at midrapidity. Two different initial conditions based 

228 on , (i) Glauber model and (ii) Color Glass Condensate (CGC), are used to model the initial 

229 transverse energy density distribution in Pb-Pb collision. The simulations are carried out for 

230 tj/s values between 0.0 to 0.16, using a lattice + hadron resonance gas model based equation 

231 of state which has a cross over temperature for the quark-hadron transition at 175 MeV. The 

232 shear viscous corrections are considered both in the evolution equations and in the freeze-out 

233 distribution function. 

234 We have compared our simulated results to the experimental data at midrapidity on the cen- 

235 trality dependence of dN c i,/dri, v'2, and V4 as a function of pj of charged hadrons measured in 

236 Pb-Pb collisions at y^NN = 2.76 TeV. The comparisons have been done for the collision central- 

237 ity classes of 10-20%, 20-30%, 30-40%, and 40-50%. We find there is nice agreement between 

238 the experimental measured dN c hjdr\ and the simulated results for both Glauber and CGC initial 

239 conditions. The V2(pr) experimental data requires a slightly lower value of tj/s for simulations 

240 with Glauber model initialization compared to the CGC based initial conditions. For a Glauber 

241 based initial condition the data prefers a tj/s value between 0.0 to 0.12, while for CGC based 

242 initial condition it prefers a value between 0.08 to 0. 16. The experimental data on v^pj) for the 

243 higher collision centralities are in general higher compared to results from ideal hydrodynamic 

244 simulations. For the more peripheral collisions of 30-40% and 40-50% centrality, the compari- 

245 son of the simulated results to the experimental data shows that the preferred 77/ s lies between 

246 0.0 and 0.08 for both Glauber and CGC based initial conditions. 

247 The current work uses smooth initial conditions. A more realistic approach would be to use 

248 a fluctuating initial condition and carry out event-by-event hydrodynamics. This will enable us 

249 to study the odd flow harmonics V3 along with the even harmonics v'2 and v\. The simultaneous 

250 description of all these experimentally measured flow harmonics in a viscous hydrodynamics 

251 framework will probably provide a better estimation of r\js at the LHC energies. We plan to 

252 consider some of these effects in the near future. The current results however agrees well with 

253 other estimates of 77/ s for the same experimental data. These estimates are based on both trans- 

254 port and hydrodynamic approaches. All results seems to indicate that the value of tj/s for the 

255 QCD matter formed Pb-Pb collisions at t/sJ^ = 2.76 TeV is between 1 - 4 x 1/47T (the KSS 

256 limit). This value is very close to that obtained for QCD matter at RHIC energies ( ysjviv = 200 

257 Gev) JHEiiiBH. 
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